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PACS ll.10.Kk - Field theories in dimensions other than four 

PACS 11 . 30 . Qc - Phase transition 

PACS 1 1 . 10 . Wx - Finite-temperature field theory 

Abstract. - We consider the 3-dimensional massive Gross-Neveu model at finite temperature as an 
effective theory for strong interactions. Using the Matsubara imaginary time formalism, we derive 
a closed form for the renormalized T-dependent four-point function. This gives a singularity, 
suggesting a phase transition. Considering the free energy we obtain the T-dependent mass, 
which goes to zero for some temperature. These results lead us to the conclusion that there is a 
second-order phase transition. 



Introduction. — The Gross-Neveu (GN) model [I], 
proposed as an effective model for quantum chromody- 
namics (QCD), has been analyzed extensively in recent 
years. In several cases, these studies consider the tempera- 
ture effect, associated with fermion systems, in both parti- 
cle and condensed matter physics [2HT8] . The motivations 
and the nature of such analysis are multiple, including in- 
vestigations of the continuous and discrete chiral symme- 
try [6l[7]. The version with N massless fermions in (2+1)- 
dimensions, for instance, presents chirality-breaking in 
perturbative analysis, and its restoration at finite temper- 
ature [14]. This provides insight into the intricate struc- 
ture of the hadronic matter, such as for the quark con- 
fmement/deconfinement transition [51 [T5HT7] . Due to such 
characteristics, the GN model is taken as a prototype for 
analysis of phase transitions in quantum field theory for 
fermions. In this sense, the search for analytical results is 
another aspect that attracts interest, in general [18], but in 
particular for N = 1 , where lattice calculations are hard to 
carry out. In this note, we address the (2+l)-dimensional 
GN model at finite temperature, but with non-zero mass 
and one type of fermion. 

The massive GN model in (2+1) dimensions was con- 
sidered in [19) . even though this version for dimensions 
greater than two is perturbatively non-renormalizable. 
This is justified by considering that, perturbative renor- 



malizability is not an absolute criterion for an effective 
theory to be physically consistent. This fact is known 
and has been used for long [2Hl25j . Here we will con- 
sider the (2+l)-dimensional massive GN model at finite 
temperature with an arbitrary zero-temperature coupling 
constant, as an effective theory for QCD. We will take 
both the thermal mass and the thermal coupling constant 
at one-loop order. This approximation will be justified 
by a non-perturbative analysis of the four-point function, 
by summing the chains of loop diagrams. This provides a 
non-perturbative relationship between the four-point func- 
tion and temperature. Such a relation leads to a singular- 
ity, suggesting a possible phase transition |26j . In order to 
get the nature of the phase transition, an analysis of the 
free energy is carried out. It gives rise to a T-dependent 
mass that goes to zero for some values of temperature, sug- 
gesting a second-order phase transition. An aditional ar- 
gument in favour of a second-order phase transition comes 
from an analysis of the beta-function. The argument relies 
on the fact that this analysis gives the result that there is 
a non-trivial infrared stable fixed point. 

Such results bring up an interesting question: can such 
a study of an effective model yield any information about a 
phase transition in physically relevant situations for QCD? 
In order to get an answer to this question, we rely on 
a procedure that follows methods first developed for the 
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BCS model in superconductivity [3S] . 

The paper is organized in the following way. In Section 
2, the model is defined. In Section 3, the thermal cor- 
rection of the mass at the one-loop level is calculated. In 
Section 4, using the temperature-dependent coupling con- 
stant, an expression for the thermal mass is derived. In 
Section 5, starting with non-perturbative results for the 
four-point function, the analysis of the critical region is 
discussed. The concluding remarks are presented in Sec- 
tion 6. 

The model. — In a D-dimensional Euclidian man- 
ifold, IR 13 , we consider the Hamiltonian for the massive 
GN model, 

H = [ d 3 x{^(x)('y j (id :j ))i(j(x)-mo^(x)i(j(x) 



(1) 



At the one-loop level the finite-temperature self-energy 
is given by 
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In order to use the dimensional regularization procedure, 
we introduce dimensionless quantities, 

Qj = kj/2TTmo , j = 1, 2, D — 1 and a — (mofi)~ 2 . (5) 

For the present case, we have D = 3 and s = 1. After 
dimensional regularization [27], we obtain 

m,T,s)= , A J%^M_, a £ (o£ + (6) 



(47r)( £, - 1 )/ 2 r(s)/3 
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where v = s — (D — l)/2. The sum in Eq. ^ is cast in 
the general form 



where mo and Ao are respectively the physical zero- 
temperature mass and coupling constant, x G M. D and 
the 7-matrices are elements of the Clifford algebra (nat- 
ural units H = c = ks = 1 are used.). This Hamiltonian 
is obtained using conventions for Euclidian field theories 
in [37]. 

From Eq. ((TJ), introducing the thermally corrected mass, 



m(T) = mo + S(T), 



(2) 



we have a free energy density of the Ginzburg-Landau 
type [33], 

T = a + b{T)<p 2 {x) +c0 4 (x) (3) 

where b{T) = —m(T) and c = Ao/2. The minus sign for 
the mass in Eq. ([3]) implies that, in the disordered phase 
we have m{T) < and for the ordered phase m(T) > 0, 
consistently. The second order phase transition occurs at 
the temperature where m(T) changes sign from negative to 
positive, characterizing a spontaneous symmetry breaking. 
In this formalism, the quantity 4>{x) = (x)ip(x) plays 
the role of the order parameter for the transition. 

Thermal self-energy. — We evaluate the thermal 
self energy, S(T), by using the Matsubara imaginary-time 
formalism. The Cartesian coordinates are specified by 
x = (x° = r, x), where x is a (D — l)-dimensional vector. 
The conjugate momentum of x is denoted by k — (fc , k), k 
being a (D — 1) -dimensional vector in momentum space. 
The KMS (Kubo, Martin, Schwinger) condition implies 
that the Feynman rules are modified by the well-known 
Matsubara prescription, [29] 
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where Lu n are Matsubara frequencies and (3 = 1 /T, T being 
the temperature. 
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is the generalized Epstein zeiet-function [30] defined for 
Re(i^) > 1/2. We then analytically continue Z\ (v,a) to 
the whole complex i/-plane 
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where K a (x) is the Bessel function of the third kind. As 
a consequence, T, s) reads 



£(AT, s) = E(D,«) 



,1-21/ 



where c = l/2w. Using Eq. 
T,(D,T,s) as 



[4»Zt c {v,a)-Z{ {u,a)\ 
(10) 

O in Eq. JTD]), we obtain 



s) = F (D, s) + Fi(D,a, s). 



(11) 



Here Fq(D, s) is a term independent of a, arising from the 
first term in brackets in Eq. ©, while Fi(D,a,s) is the 
term arising from the sum in the second term. We find 
that, for s = 1 and even dimensions D > 2, Fo(D,s) is 
divergent due to the pole of the T-function. Accordingly, 
this term is subtracted to get the physical a-dependent 
function S^(Z?,a). The mass counter-term is a pole ap- 
pearing at the physical value s = 1. The a-dependent 
correction to the mass is proportional to the regular part 
of the analytical extension of the Epstein zeia-function in 
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the neighborhood of the pole at s = 1 . For uniformity, the 
term F (D) is also subtracted in the case of odd dimen- 
sions D, where no poles of T-functions are present. There- 
fore, using the modified Bessel function for D = 3, s = 1 



and a = (m /T)~ 
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physical thermal self-energy for D = 3 is given by 
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The behavior of f(T) is the following: for T — > 0, /(T) — > 
0; for T — > oo, /(T) -> -oo, with /(T) < for all values 
of T > 0. Using Eq. the temperature dependent mass 
is 

(14) 



m(T) = TOo + _A m2/(T). 

The condition for a phase transition, m(T) — 0, provides 
a critical temperature for each fixed value of the zero- 
temperature coupling constant, Ao- This result indicates 
a second-order phase transition. However, in order to es- 
tablish this result on a firmer ground, we consider the 
T-dependent correction to the coupling constant Ao- 

T-dependent coupling constant. Initially, we 
consider the four-point function with null external mo- 
menta, which defines the temperature-dependent coupling 
constant. It is given up to one- loop by, 



r^ 4) (A ,/3)~ A [l + A n(A/3)] , 



(15) 



where U(D,f3) is the /3-dependent one-loop polarization 
diagram given by 
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Using the dimensionless quantities introduced in Eq. ([5]), 
H(D, 0) is given as 
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The function U d a) is extended to the whole complex 
fi -plane, resulting in 
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where h(fi,D) = n 2 ^- D / 2 /(2 D - 2 T(fi)) and 
W(p,a,D) = 2^ 



n=l 



Using these results in Eq. (|17p . we obtain II(£>, a, s) as 

II(D,a,s) = +G(D,a,s), (21) 

where H(D, s) is a term independent of a, arising from the 
first term in brackets in Eq. (fT9|) . while the term G(Z?, a, s) 
arises from the the TU-functions. As in the case of the 
thermal self-energy, the physical thermal polarization is 
obtained by subtracting the term H(D, s) in Eq. (|2"Tj) and 
is given by, 
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Fig. 1: Polarization (in units of mo), S = I1h(3, T)/mo, as a 
function of T (in units of mo, i.e. t — T /mo). 



(17) Taking D = 3 and a = (ra /T)~ 2 , we obtain 



n*(3,T) = ^ 



— In ( 1 

m 



-mo/T 



I 



1 + e m /T 



. (23) 

Substituting this expression into Eq. (I15|) . we obtain the 
temperature-dependent coupling constant, 

g(T; A ) = r^(A ,T) ~ A [l + A n R (3,T)]. (24) 

The physical polarization Hr (3, T) is shown in Fig. [1] We 
find from this plot and from Eq. flM|) that g(T = 0; A ) = 
Ao; this just reflects the fact that Ao is the physical cou- 
pling constant at T = 0. 
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The temperature-dependent coupling 
g(T;Xo), is introduced in place of Ao in 
Then we get a T-dependent mass, 

m(T) - m [l + ^-g(T;\ )mof(T)]. 

Z7T 



constant, 
Eq. (HI. 



(25) 



We have then all the elements to analyze the phase tran- 
sition. 

The phase transition. — Let us first remind an es- 
sential feature of a situation in condensed matter physics, 
theoretically close to the one we examine in this article; 
namely, the BCS field theoretical approach to the super- 
conducting phase transition |26j . In this case it is shown 
that, at criticality, the leading contribution to the four- 
point function with zero-external momenta is given by 
the sum of all chains of one-loop diagrams. This non- 
perturbative calculation leads to an expression of the form 



[l-A n fl (3,T)] 



(26) 



In our case the first two terms of the expansion in powers 
of Ao of such a function are given in Eq. (|2"4"|) . The exis- 
tence of a singularity of the four-point function in Eq. ([26| 
indicates a phase transition, as explained in |26j ; in other 
words the phase transition is characterized by the poles of 
ifj(A ,T c ),i.e., 



A IM3,T) = 1. 



(27) 



In our case this singularity is found from an analysis of 
Fig. [1] The nature of the transition is obtained by a study 
of the free energy, Eq. ([3]). 

In the present article we have a situation analogous to 
the BCS theory, and we have done the study of the free 
energy for a second-order phase transition starting from 
Eq. ©. Then the the critical region is defined by taking 
m(T c ) = 0, in Eq. ([25]), which leads to 



l(T c ; A ) 



2n 



m f(T c ) ' 



(28) 



This provides T c as a function of Ao, giving a curve plot- 
ted in Fig. [5] On the curve we have m(T c ) = 0; below the 
curve we have m(T) > and above the curve m(T) < 0. 
We find that the sign for the thermal mass corresponds 
to the expected behavior of a phase transition from a dis- 
ordered to an ordered phase, as temperature is lowered. 

Moreover, Eq. (|2"T)) gives rise to a function t — t(X) 
(where again A = Aomo/27r and t — T/mo), with a be- 
havior that is similar to that presented in Fig. [2j Indeed, 
from Fig. Q] we obtain that t — > 0.87 for Ao — > oo. On the 
other hand, t — ¥ oo for A — > 0. This non-perturbative re- 
sult is then compatible with the analysis performed in the 
previous sections, of the mass as a function of the temper- 
ature, taking the thermal corections to both the mass and 
the coupling constant at one-loop level. Then we conclude 
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Fig. 2: Phase transition defined by m(T c ) = 0: values for 
the pair (T c ,Ao) define the critical line; on the vertical axis 
t = T c /mo and on the horizontal axis A = Aomo/27r. The 
minimum critical temperature is T c ~ 0.874mo. 



that indeed there is a second-order phase transition, char- 
acterized by the divergence of the four point function and 
with the critical temperature obtained from the condition 
m{T c ) = 0. 

A further indication of a second-order transition can be 
obtained from a renormalization group point of view, fol- 
lowing lines analogous to those employed in |35j for type-II 
superconducting films. If an infrared stable fixed point ex- 
ists, it is possible to determine it by a study of the infrared 
behaviour of the beta-function, i.e, in the neighbourhood 
of vanishing external momentum \p\ = 0; we consider the 
thermal coupling constant at criticality, with an external 
small momentum, given by 



g(T, \p\ « 0) 



Ao 



[1-A IIOD,7>)] 



(29) 



Using a Feynman parameter x, the one-loop four-point 
function is, 
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The second term in Eq. (f30|) vanishes in the limit \p\ — > 
0. Following the calculation steps for small values of \p\ 
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detailed in [35]; we are left for \p\ ps 0, after regularization, 
with an expression similar to Eq. (j2"2")l . 
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U R (D,T,p) = 2D _° 2nD/2 [2nW(2,a c (p),D) 
-W(l,a c (p),D)}, 



(32) 



where it is understood that we take asymptotic values for 
small values of the argument in the Bessel functions in 
Eq. (f20|) . which defines the quantity W(fi,a,D) . 
For D = 3 we have, M®~ 2 ps to + \p 2 x(l - x), and 



Il R (3,T,p)=B(3,T)+A{3,T)\ P \ 



(33) 



where, 



to 



B ( 3 ' T ) = ^2 l 2 *W(2, a c , 3) - W(l, oc, 3)] (34) 



and 
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2tt 3 / 2 
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with a = Jq 1 dccy / a;(l — x) = n/8. 

The coupling constant has dimension of \p\~ , taking 
\p\ as a running scale we define a dimensionless coupling 
constant 



g' = \p\g = 



1-X [B(T, 3]+A(T,3)|p| 



and the beta-function, 



(36) 



we easily see that the condition of a non-trivial infrared 
stable fixed point is fullfilcd by the solution 
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A(3,T C 



(37) 



As we have already stated in this note, the GN model 
may be seen as an effective theory for QCD, with an arbi- 
trary zero-temperature coupling constant. Let us then es- 
timate a specific value for the critical temperature, choos- 
ing the mass of the Gross-Neveu fermion to be the effective 
quark mass of the proton [36] , mo ps 330 MeV. Consider a 
strong coupling regime characterized by large values of Ao , 
up to the limit Ao — > oo (see Fig. [5]). For these values (for 
instance for Ao ~ 167t/too = 0.15 MeV -1 , corresponding 
to A = 8 in Fig. [2]), we find a critical temperature, T c , of 
the order of T c ~ 0.87to ps 288 MeV. 

Concluding remarks. — In short, we have answered 
positively the question whether there is a phase transi- 
tion in the N = 1 massive GN model; i.e. yes, there is 
a second-order phase transition. This conclusion is based 
in two results. From the non-perturbative analysis of the 
four-point function, a critical region can be established. 



Thereby the analysis of the free-energy provides the na- 
ture of the transition, as being of second-order. It is im- 
portant to add that this procedure is strikingly similar 
in character to the method for finding critical tempera- 
ture for the second-order phase transition in BCS model 
of superconductivity. Also, we have shown that we can de- 
fine a freta-function, which has a non-trivial infrared fixed 
point, thus reinforcing the conclusion that the transition 
is a second-order one. 

Considering a regime of strong coupling as defined 
above, and the mass of the GN fermion as the effective 
quark mass, we find a critical temperature of the order of 
T c ~ 288MeV. It is interesting to observe that this value 
is of the same order of magnitude as the estimated temper- 
ature ~ 200MeV for the quark deconfinement transition, 
obtained from lattice calculations. The phase-transition 
here may be then associated with the transition from a 
hadronic state to a quark-gluon plasma. 

This has been possible by employing the massive GN 
model in three dimensions, which has been shown to ex- 
ist and constructed for the first time in |19j . In previous 
works devoted to get insights on the behaviour of hadronic 
matter [2 13, T5"HT7] . the massless GN model, in its version 
with a large number of components, has been often em- 
ployed. Using the one-component massive GN model in 
three dimensions and taking the fermion mass as a phys- 
ical parameter (the effective quark mass), we have been 
able using analytical means, to determine the transition 
temperature. Moreover this temperature is found to fall 
into a range of values compatible with the transition tem- 
perature for hadronization, estimated from lattice simula- 
tions. A rigorous study involving aditional aspects of this 
transition, including determination of critical exponents, 
is left for future work. 

Finally, we would like to add some comments about the 
suitability of the GN model to fit QCD. We have in mind 
the Nambu-Jona-Lasinio (NJL) model which, particularly 
in its four-dimensional version is closer to QCD than the 
GN model. Many studies on the phase structure of the two 
and four-dimensional NJL model have been performed, 
which includes articles by some of us (see [371[35] and ref- 
erences therein). A relevant aspect of the massive three 
dimensional GN model, as we have employed in this arti- 
cle, is that the critical temperature can be obtained ana- 
lytically in the context of a well known phenomenological 
approach to phase transitions. This is much harder to be 
studied with the NJL model; it will the subject of future 
work. 
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